###
### NTHROOT.R  +++ Test suite +++
###


test.nthroot <- function(input, expected) {
    output <- do.call(getFromNamespace("nthroot", "matlab4r"), input)
    identical(output, expected)
}

nthroot.expected.n0 <- c(-1.2)
nthroot.expected.n1 <- c(1, -2, 3)
nthroot.expected.n2 <- c(1,  0, 3)
nthroot.expected.n3 <- c(1, -2, 3)

test.nthroot(list(x=-1.2^5, n=5), nthroot.expected.n0)
test.nthroot(list(x=c(1,-2, 3), n=1), nthroot.expected.n1)
test.nthroot(list(x=c(1, 0, 9), n=2), nthroot.expected.n2)
test.nthroot(list(x=c(1,-8,27), n=3), nthroot.expected.n3)
